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Field strength correlators in QCD at zero and non-zero temperature 
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We study, by numerical simulations on a lattice, the behaviour of the gauge-invariant field strength correlators 
in QCD both at zero temperature, down to a distance of 0.1 fm, and at finite temperature, across the deconfinement 
phase transition. 



1. INTRODUCTION 

An important role in hadron physics is played 
by the gauge-invariant two-point correlators of 
the field strengths in the QCD vacuum. They are 
defined as 



r (x) = (0\Tr{G„ p (x)SG ua (0)S^} |0) 



(1) 



where G MP = <7T a G° p is the field-stength tensor 
and S = S(x,0) is the Schwinger phase operator 
needed to parallel-transport the tensor G ua (0) to 
the point x. 

These correlators govern the effect of the gluon 
condensate on the level splittings in the spectrum 
of heavy QQ bound states [[j] || . They are the ba- 
sic quantities in models of stochastic confinement 
of colour |^-^| and in the description of high- 
energy hadron scattering J7| — [3~0t] - 

A numerical determination of the correlators 
on lattice (with gauge group ££7(3)) already ex- 
ists, in the range of physical distances between 



0.4 and 1 fm |lj 
exponentially 



In that range V 



falls off 



r (x) ~ exp(-\x\/X) 



with a correlation length A ~ 0.22 fm | Ejj . 

What makes the determination of the corre- 
lators possible on the lattice, with a reasonable 
computing power, is the idea (l^Jl^l of remov- 
ing the effects of short-range fluctuations on large 
distance correlators by a local cooling procedure. 

* Speaker at the conference. 



Freezing the links of QCD configurations one after 
the other, damps very rapidly the modes of short 
wavelength, but requires a number n of cooling 
steps proportional to the square of the distance d 
in lattice units to affect modes of wavelength d: 



kd 2 



(3) 



Cooling is a kind of diffusion process. If d is suf- 
ficiently large, there will be a range of values of 
n in which lattice artefacts due to short-range 
fluctuations have been removed, without touching 
the physics at distance d; by lattice artefacts we 
mean statistical fluctuations and renormalization 
effects from lattice to continuum. This removal 
will show up as a plateau in the dependence of 
the correlators on n. This was the technique suc- 
cessfully used in Ref. pd] |. There, the range of 
distances explored was from from 3-4 up to 7-8 
lattice spacings at (3 ~ 6 (/3 = Q/g 2 ), which means 
approximately from 0.4 up to 1 fm in physical dis- 
tance. The lattice size was 16 4 . 

In Sect. 2 we discuss new results obtained on 
a 32 4 lattice, at (3 between 6.6 and 7.2: at these 
values of j3 the lattice size is still bigger than 1 fm, 
and therefore safe from infrared artefacts, but d = 
3,4 lattice spacings now correspond to physical 
distances of about 0.1 fm. Since what matters 
to our cooling procedure is the distance in lattice 
units, we obtain in this way a determination of 
the correlators at distances down to 0.1 fm JlJ]. 

All of the above concerns the theory at zero 
temperature. In Sect. 3 we go further and de- 



2 



termine the behaviour of the correlators at fi- 
nite temperature for the pure-gauge theory with 
SU(3) colour group and in particular we study 
their behaviour across the deconfining phase tran- 
sition pa . The motivations to do that stem from 
Refs. fflil. 



2. RESULTS AT T = 



The most general form of the correlator com- 
patible with the invariances of the system at zero 
temperature is 



~\~\Xfj i Xi/Sp(j x px^Spi, -\~ x pX(j8 py 



XpXy Spa- ) - 



dx 2 



(4) 



T> and T>\ are invariant functions of x 2 . We work 
in the Euclidean region. 

It is convenient to define a X>||(a; 2 ) and a 
D±(x 2 ) as follows: 



V + V l + x 

v + v x . 



,&Di 
dx 2 



(5) 



On the lattice we can define a lattice operator 
Trt p va , which is proportional to T> WtVcr in the 
naive continuum limit, i.e., when the lattice spac- 
ing a — > 0. Making use of the definition (5) we 
can thus write, in the same limit, 



V\{da) ~ o a 4 £>||(dV) + 0(a 6 ) . 
2>£(da) ~ a A V ± {d 2 a 2 ) + 0{a 6 ) 

a— »Q 



(6) 



Equations (6) come from a formal expansion of 
the operator, and are expected to be modified, 
when the expectation value is computed, by lat- 
tice artefacts, i.e., by effects due to the ultraviolet 
cutoff. These effects can be estimated in pertur- 
bation theory and subtracted . Instead we re- 
move them by cooling the quantum fluctuations 
at the scale of the lattice spacing, as explained in 
the Introduction. After cooling, T>^ and T>\ are 
expected to obey Eqs. (6). The typical behaviour 
of T>^ and of T>\ along cooling is shown in Fig. 
1. Our data are the values of the correlations at 
the plateau. 



COOLING STEP 



Figure 1. A typical behaviour of T>^ (diamonds; 

d = 6, = 6.6, lattice 32 4 ) and of T>\ (squares; 
d = 12, (3 = 6.6, lattice 32 4 ) during cooling. 



We have measured the correlations on a 32 4 
lattice at distances ranging from 3 to 14 lattice 
spacings and at (3 = 6.6, 6.8, 7.0, 7.2. From 
renormalization group arguments, 



a((3) 



1 

a7 



/(/?) , 



(7) 



where A^ is the fundamental mass-scale of QCD 
in the lattice renormalization scheme. At large 
enough (3, f((3) is given by the usual two-loop 
expression: 



51/121 



exp 



33 M 



(8) 



for gauge group 5/7(3) and in the absence of 
quarks. At sufficiently large (3 one also expects 
that 



2t/(/3)- 4 
2?i/(/3)- 4 



1 



(9) 



where f{(3) is given by Eq. (8) and terms of higher 
order in a are negligible. 

In Fig. 2 we plot £>|/(/3)- 4 and V\f{{3)-^ ver- 
sus dphys = (d/h-if) /(/?)• I n this figure we have 
also plotted the values of the correlators obtained 
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Figure 2. The functions T>j_f(/3) 4 (upper curve) 
and Djj"/(/3)~ 4 (lower curve) versus physical dis- 
tance (in fermi units). Triangles correspond to 
the data of Ref. [11]. The lines are the curves 
for T>± and T>u obtained from the best fit of Eqs. 
(10) and (11). 



in Ref. jf]J, corresponding to physical distances 
rfphys > 0.4 fm. We have applied a best fit to all 
of these data with the functions 



V (1) {x 2 ) = A w exp(-\x\/X A ) 
+%jl ex p(-\x\/\ a ) . 

We have obtained the following results: 



(10) 



A 



3.3 x 10 s 

a ~ 0.69 
1 1 
~ 7U 182 



-j- ~ 0.7 x 10 s 

ai ~ 0.46 , 
1 1 



A. 



(11) 



with x /-Wd.o.f. — 1-7- The continuum lines in 
Fig. 2 have been obtained using the parameters 
of this best fit. With the value of Al determined 
from the string tension |B(I we obtain 



X A ~ 0.22 fm , A a ~ 0.43 fm 



(12) 



The correlation length Xa, which enters the non- 
perturbative exponential terms of T> and X>i, as 
well as the magnitude of the coefficients A and 
A±, are compatible with the values obtained in 
Ref. Jill. 



We stress again that we have been able to ob- 
serve terms proportional to l/|a;| 4 in the corre- 
lations because we have worked at larger values 
of /3, where the distance between two points (far 
enough in lattice units so that the correlation is 
not modified by cooling before lattice artefacts 
are eliminated) is small compared with 1 fm in 
physical units. A larger lattice (32 4 ) has been 
necessary to avoid infrared artefacts. 

3. RESULTS AT FINITE T 

To simulate the system at finite temperature, a 
lattice is used of spatial extent N a » N T , N T be- 
ing the temporal extent, with periodic boundary 
conditions. The temperature T corresponding to 
a given value of (3 = 6/g 2 is given by 

N T ■ a((3) = i , (13) 

where a{0) is the lattice spacing, whose expres- 
sion in terms of (3 is given by Eqs. (7) and (8). 

At finite temperature, the 0(4) space-time 
symmetry is broken down to the spatial 0(3) 
symmetry and in principle the bilocal correlators 
(1) are now expressed in terms of five indepen- 
dent functions [p~6| p~8|| (instead of two as in the 
zero-temperature case); two of them are needed 
to describe the electric-electric correlations: 

(0\Tr {E i (x)S(x,y)E k (y)SHx,y)} |0) = 
,dDf 



Oik 



4 du\ 



dDf . N 



where = is the electric field operator and 
u-fi = - [u 2 = (x- y) 2 }. 

Two further functions are needed for the 
magnetic-magnetic correlations: 

(0\Tr {B i (x)S(x,y)B k (y)S^x,y)} |0) = 
t dD^ 



Hk 



D B +Df+ u 2 



dip 



dDf . . 

UlUk ~m 1 (15) 



where B k = ^eijkGij is the magnetic field oper- 
ator. Finally, one more function, D BE , is neces- 
sary to describe the mixed electric-magnetic cor- 
relations. 

The five quantities D E ', Df , D B , Df and D BE 
are all functions of u 2 , due to rotational invari- 
ance, and of u 2 , due to time-reversal invariance. 



4 



From the conclusions of Refs. (T^-[l8), one ex- 
pects that D E is related to the (temporal) string 
tension and should have a drop just above the 
deconfinement critical temperature T c . In other 
words, D E is expected to be the order parameter 
of the confinement. Similarly, D B is related to 
the spatial string tension |l6| , [l7|| . 

We have determined the following four quanti- 
ties 



Vf(u 2 ,0) 
V E (u 2 ,0) 
V B (u 2 ,0) 
V B (u 2 ,0) 



V E (u 2 ,0) +V E (u 2 ,0) +u 
V E (u 2 ,0)+V?(u 2 ,0) 
V B (u 2 ,0) +Vf (u 2 ,0) +u 
V B (u 2 ,0)+V B (u 2 ,0) , 



dip 



• d-D? 
dip 
(16) 




y 1.05 T/Tc 



Figure 3. The quantity Vf/A\ [defined by the 
first Eq. (16)] versus T/T c and the physical spa- 
tial distance (in fm). 



by measuring appropriate linear superpositions 
of the correlators (14) and (15) at equal times 
(114 — 0), on a 16 3 x4 lattice (N T = 4, in our nota- 
tion). The critical temperature T c for such a lat- 
tice corresponds to (3 C ~ 5.69. The behaviour of 
T> E and T> E is shown in Figs. 3 and 4 respectively, 
on three-dimensional plots versus T/T c and the 
physical spatial distance. Due to the logarithmic 
scale, the errors are comparable with the size of 
the symbols and the lines connecting the points 
are drawn as an eye-guide. A clear drop is ob- 
served for T> E and V E across the phase transition, 
as expected. 

On the contrary, no drop is visible across the 
transition for the magnetic correlations T> B and 
T> B : as an example, the behaviour for T> B is 
shown in Fig. 5. 

Our results can be summarized as follows: 

(i) In the confined phase (T < T c ), until very 
near to the temperature of deconfinement, the 
correlators, both the electric-electric type (14) 
and the magnetic-magnetic type (15), are nearly 
equal to the correlators at zero temperature (m) : 
in other words, D E ~ D B ~ D and Df ~ Df ~ 
D 1 for T <T C . 

(ii) Immediately above T c , the electric-electric 
correlators (14) have a clear drop, while the 
magnetic-magnetic correlators (15) stay almost 
unchanged, or show a slight increase. More pre- 



cisely, looking at the values for the difference 



V E (u 2 ,0) 



V E (u 2 ,0) 



BT) E 



dip 



(17) 



between the two quantities reported in Figs. 4 
and 3 respectively, one finds that the quantity T> E 
does not show any drop across the phase transi- 
tion at T c . So the clear drop seen in the quantities 
T> E and V E across T c is entirely due to a drop of 

D E alone. This result confirms the conclusion of 
Refs. where V E was related to the (tem- 

poral) string tension erg. Similarly, one can look 
at the following difference: 



V B (u 2 1 Q)-V B (u 2 ,0) 



.^^(^o). (18) 



dip 



One thus finds that T>f does not show any drop 
across the transition and, in addition, it is nearly 
equal to Vf (Df ~ Vf ). From the fact that the 
quantities T> B and T> B stay almost unchanged (or 
even show a slight increase) across T c , we con- 
clude that the same must be true also for V B . 
It was shown in Refs. fl6|,[l7j that V B is related 
to the spatial string tension a s . Recent lattice 
results |2lj indicate that <r s is almost constant 
around T c and increases for T > 2T C : this fact 
is in good agreement with the behaviour that we 
have found for T> B . 
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